RF Shimming Pulses For Ex-Situ NMR Spectroscopy and 
Imaging Using Bl Inhomogeneities 
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I describe a method for generating "shim pulses" for NMR spectroscopy and imaging (MRI) by 
taking advantage of the inherent inhomogeneity in the static and radiofrequency (RF) fields of a 
one-sided NMR system. The RF inhomogeneity here is assumed, without loss of generality, to be a 
linear gradient. General polynomials in the spatial variables can be generated using x, y and z RF 
gradients using trains of hard pulses which result in linear combinations of monomials xy, y 2 , xz 
etc., and any desired scalings of these monomials. The basic shim pulse is constructed using small 
tip angle approximations. 

PACS numbers: 76.60.Jx 



I. INTRODUCTION 

"Shim pulses" have been proposed by Topgaard and 
Pines [l[ using pairs of adiabatic pulses in combination 
with gradient modulations. Shim pulses for spectroscopy 
are meant to be used intermittently during stroboscopic 
FID readout; they work by correcting the phase of spins 
as function of spatial position while leaving the chemi- 
cal shift evolution intact. Classical Hahn echo refocusing 
pulses are obviously inadequate since they would refocus 
the chemical shift as well, removing all spectral informa- 
tion. 

For imaging purposes, shim pulses can be used to im- 
part the correct phase at the center of k-space acquisi- 
tion. The main challenges to designing shim pulses are: 
1) keeping the pulse as short as possible so that large 
spectral widths can be used, 2) finding the best way to 
impart a large enough phase modulation given the lim- 
ited available gradient strength during this short amount 
of time. 

The adiabatic pulses of Topgaard are generally too long 
for many applications and this limits the performance 
of many experiments. For example, a 10 ms delay be- 
tween consecutively collected samples results in a maxi- 
mum spectral width of 100 Hz, and readouts of several 
points lead to significant T2 relaxation. 

Here, we exploit the idea that in an ex-situ NMR en- 
vironment (e.g. with a single-sided NMR magnet con- 
figuration) , the static field is intrinsically inhomogeneous 
and so is the radiofrequency (RF) field due to the nature 
of single-sided magnet and coil designs. It thus makes 
sense to think of schemes that take advantage of the in- 
herent RF inhomogeneity, available at no additional cost, 
to correct for the negative effects of static field inhomo- 
geneities. 

I present a method which uses trains of hard pulses 
generated by inhomogeneous RF fields. It is based on 
the approach of Meriles et al. Q in the sense that hard 
pulses are used. However, we explicitly construct the de- 
sired polynomials. In the limit of small flip angles, these 
rotations commute and combine to create any desired 



polynomial in the spatial variables x, y and z. The cal- 
culations presented herein demonstrate the ability of RF 
gradients to generate z-rotation shim pulses or excitation 
pulses. 



II. THEORY 



Consider the following product of four rotations: 



e x e Y e 



■X„-Y 



(1) 



where the dots mean "higher order terms" . In what fol- 
lows, we drop these higher order terms from the notation. 

We observe that if X contains the spin operator I x 
and Y contains I y , the commutator [X, Y] will contain 
I z . This principle allows us to generate a rotation about 
the z axis. Moreover, the product of X and Y in this 
commutator allows us to generate polynomials in the co- 
efficients of X and Y. 

We consider therefore rotation operators of the form: 



Xi = giXiI x , Yj = hjijjly (2) 

where Xi and j/j refer to any of the spatial variables x, y 
or z. The commutator of Xj and Yj is: 

[Xj,Y j }=ig j h j x j y j I z (3) 

and allows us to generate any of the following gradient 
terms in the monomials x 2 , y 2 , z 2 , xy, yz or xz by pick- 
ing Xj and yj appropriately. Moreover, we may pick the 
product gjhj to produce any desired scaling (rotation an- 
gle). 



A. Second order shim pulse 

We may combine N successive groups of 4 pulses back 
to back to generate any linear combinations of N mono- 
mials as follows: 
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We will call this type of pulse a second order shim pulse 
(because it contains the product Xj-yj of spatial variables) 
and denote the basic unit e X:i e Yj e~ Xj e~ Yj as S 2 {Xj, Yj). 

While this 5*2 pulse creates an I z rotation, it is ideal 
for use in stroboscopic pulse train experiments, where the 
phase of freely evolving spins is periodically corrected. 
The S2 pulse can be converted into an excitation pulse 
using the following method: the 5*2 pulse is sandwiched 
on the left by a (iv/2)y on the left and (ir/2)~ y on the 
right. This has the effect of rotating I z into I x , which can 
then serve as an excitation pulse. If multiple subunits of 
S2 are used, the sandwiching need only be applied to 
the combined group of pulses. We note that the (ir/2) 
pulse should be applied using a homogeneous RF field or 
using composite pulses that compensate for the spatial 
inhomogeneity of the RF field. 

B. Third order shim pulse 

A third order shim pulse can be generated from the 
following sequence: 

S 2 (X, Y)e z S 2 (X,Y)e- z = e ^ Y ^ = S 3 (X, Y, Z) (5) 

This double commutator is unable to generate an I z 
rotation if the building blocks are the operators I x and 
I y . In this case, the S3 pulse subunit generates an I x or 
I y rotation which can serve as excitation pulse. 

To get an I z rotation pulse, the following modification 
to S3 can be used: an I x rotation can be converted into 
an I z rotation using a 7r/2 rotation about I y . There- 
fore, sandwiching the S3 pulse, by (tt/2) v on the left and 
(n/2)~ y on the right converts it to an I z rotation. We 
note that the (7r/2) pulse should be applied using a ho- 
mogeneous RF field or using composite pulses that com- 
pensate for the spatial inhomogeneity of the RF field. 

C. Fourth order shim pulse 

Generalizing these ideas, we may create a fourth order 
shim pulse can using the following sequence: 

S 3 (X,Y,Z)e w S 3 (X,Y,Z)e- w = 

e w,r\,z\w = Si (x,Y,z,w) (6) 

This sequence of commutators can generate an I z ro- 
tation, for example, as follows: 

[[[I x ,I y ],I x ],I x ] = -iI z (7) 



D. Higher order shim pulses 

A fifth order shim pulse is given by: 



Si(X, Y, Z, W)e u X 1 Y, Z, We~ u = 

e [[[[X,Y],Z},W},U] 

= S 5 (X,Y,Z,W,U) (8) 
whereas a sixth order shim pulse is the sequence: 



S 5 (X, Y, Z, W, U)e v S 5 {X, Y, Z, W, U)e~ v = 

e [[[[[X,Y],Z],W],U],V] 

= S 6 (X,Y,Z,W,U,V) (9) 

This latter pulse generates an I z rotation, for example, 
as follows: 

[[[[[I x ,I y },I x ],I x ],I x ],I x ] = -iI z (10) 

This process of can be generalized to any nth order 
monomial of the form XjUj where a + b = n. By concate- 
nating any number of monomial units, we may generate 
any linear combinations of them, by choosing the RF 
pulse amplitudes to match the desired scalings of each 
monomial. This results in an arbitrary polynomial shim 
pulse. 

We note that the higher the order of the mononomial 
term, the higher the error in the first order approxima- 
tion to the rotation will be. For example, a monomial 
such as x 2 grows rapidly with increasing distance x from 
the origin. Therefore, the errors grow larger close the 
sample edges. One possible solution to this problem is to 
repeat the S2 sequence many times for smaller flip angles. 
The error is then reduced to any desired order provided 
the flip angles are made small enough. In the following 
section, I present results that illustrate this refinement 
approach. 

A second possible approach would be to take a basic 
sequence such as S2, repeat it a few times, and use op- 
timal control algorithms (such as GRAPE) to find the 
proper spatial variables, spin operators and their relative 
phases in order to generate a propagator with far fewer 
errors in the spatial and spin operator orders. The free 
parameters which could be used in an optimal control 
optimization include pulse phase, gradient direction (x, 
y or z) and pulse amplitude. Another method would be 
to use the inverse scattering transform to solve for the 
RF waveform given the desired excitation profile. These 
methods will not be discussed here. 



III. THE REFINEMENT METHOD 

Figure[T]shows simulations of a second order shim pulse 
x 2 + y 2 . Figure [T^a) shows the ideal or "target" phase 
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(a) desired phase (b) single unit 




FIG. 1: Simulation of a second order x 2 + y 2 shim pulse: (a) 
the desired or "target" phase profile, (b) attempt at the shim 
pulse using a single unit of S2{X,Y), (c) five units and (d) 
ten units. 

profile in the xy plane. The pulse unit required to gen- 
erate this kind of shim pulse is 

S 2 {X l ,Y 1 )S 2 (X 1 , Yi ) =e xI * e xI y e~ xI * e~ xI y 

This basic unit is the N = 2 case. It requires a total of 
8 hard pulses. The phase obtained from a single unit of 
this shim pulse is shown in Figure [Ijb) . We see that the 
pulse performs well at the origin but suffers from extreme 
distortions at large values of x 2 + y 2 . This pulse can be 
concatenated into five "smaller' (weaker amplitude) but 
otherwise identical RF pulses, as shown in Figure [ljc) 
with a slight improvement near the center. 

The ten units version of this pulse is shown in Fig- 
ure [Hd) and reproduces well the desired (ideal) phase 
pattern of Figure [TJa). This pulse contains a total of 



8 x 10 = 80 hard pulses. If each hard pulses is 10 (is 
long, the total length of this pulse would be 800 (is. A 
800 (is shim pulse is already an improvement by an order 
of magnitude over the Topgaard adiabatic pulses. 

In practice, the total duration of each hard pulse de- 
pends on the available maximum B\ amplitude. When 
considering the time-optimality of this shim pulse, only 
the pulse area for each hard pulse matters (duration x 
amplitude). Thus, in order to halve the total duration of 
the shim pulse, we require the doubling of B\. 



IV. CONCLUSION 

The basic idea of using RF gradients to generate arbi- 
trary order shim pulses works, as shown in the previous 
simulations and theoretical expressions. Two basic errors 
arise in these analytical expressions: for the generation 
of an I z shim pulse, the presence of I x and I y which "con- 
taminate" a desired I z rotation, and the actual coefficient 
of I z contains not only the desired monomial but mono- 
mials which are higher order in the spatial variables. A 
similar argument applies to I x or I y pulses for excitation 
purposes. 

It is likely that the performance of such pulses could 
be improved by optimal control methods, where the spa- 
tial variables Xj , yj , Zj pulse amplitudes gj , h j and pulse 
phases I x vs. I y (or intermediate phases) are used as op- 
timization parameters during the search. This requires 
that some fidelity measure for the propagator is used that 
minimizes the errors in the "outer parts" of the spatial 
variables, where the pulse errors are largest. By allowing 
arbitrary choices of the spatial variables and their coef- 
ficient in the optimization algorithm, it could be possi- 
ble to eliminate the higher order spatial terms (monomi- 
als) in I z via cancellations (negative terms could possibly 
cancel positive terms). By allowing the phase of the RF 
pulse elements to vary, additional desirable elements such 
as robustness to offsets, and elimination of spurious spin 
operator terms I x , I y from the propagator of a desired I z 
rotation would be possible. 
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